CMB anisotropies from vector perturbations in the bulk 
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The vector perturbations induced on the brane by gravitational waves propagating in the bulk 
are studied in a cosmological framework. Cosmic expansion arises from the brane motion in a non- 
compact five- dimensional anti-de Sitter spacetime. By solving the vector perturbation equations in 
the bulk, for generic initial conditions, we find that they give rise to growing modes on the brane in 
the Friedmann-Lemaitre era. Among these modes, we exhibit a class of normalizable perturbations, 
which are exponentially growing with respect to conformal time on the brane. The presence of 
these modes is, at least, strongly constrained by the current observations of the cosmic microwave 
background (CMB). We estimate the anisotropies they induce in the CMB, and derive quantitative 
constraints on the allowed amplitude of their primordial spectrum. Our results provide stringent 
constraints for all braneworld models with bulk inflation. 

PACS numbers: 04.50.+h, ll.10.Kk, 98.80.Cq 



I. INTRODUCTION 

The idea that our universe may have more than three 
spatial dimensions has been originally introduced by 
Nordstrom Q, Kaluza Q and Klein 0. The fact that 
super string theory, the most promising candidate for a 
theory of quantum gravity, is consistent only in ten space- 
time dimensions (11 dimensions for M-theory) has led to 
a revival of these ideas [J, la, |6| . It has also been found 
that string theories naturally predict lower dimensional 
"branes" to which fermions and gauge particles are con- 
fmed-while gravitons (and the dilaton) propagate in the 
bulk 0, IS • Such "braneworlds" have been studied in a 
phenomenological way already before the discovery that 
they are actually realized in string theory 0, |nj . 

Recently it has been emphasized that relatively large 
extra-dimensions (with typical length L ~ /jm) can 
"solve" the hierarchy problem: The effective four dimen- 
sional Newton constant given by G 4 oc G/L N can be- 
come very small even if the fundamental gravitational 
constant G ~ m^ 2+N ^ is of the order of the electro- 
weak scale. Here N denotes the number of extra- 
dimensions m m q . It has also been shown that 
extra-dimensions may even be infinite if the geometry 
contains a so-called "warp factor" [l6| . 

The size of the extra-dimensions is constrained by the 
requirement of recovering usual four-dimensional Ein- 
stein gravity on the brane, at least on scales tested by ex- 
periments [l7l llM|]ll . Models with either a small Planck 



mass in the bulk ll2Lll3lll4| , or with non-compact warped 
extra-dimensions |l5l [lfj , have until now been shown to 
lead to an acceptable cosmological phenomenology on the 
brane [S3, H3, H3, H3, HE I: 



|26j , with or without Z2 sym- 
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metry in the bulk p^ . I28L l29l | . Explicite cosmological sce- 
narii leading to a nearly Friedman-Lemaitre universe at 
late time have been proposed by means of a static brane 
in a time dependent bulk l31l or with a moving brane 
in an anti-de Sitter bulk \.'>'A |33J . It has been shown that 
both approaches are actually equivalent [34|. 

One can also describe braneworlds as topological de- 
fects in the bulk EE EE El EE E3 This approach 
is equivalent to the geometrical one in the gravity sec- 
tor |40|. while it admits an explicite mechanism to con- 
fine matter an d gaug e fields on the brane EE 0, EE 
EE E3, EH 123, |43 E3, Ull- Although depending 

on the underlying theory, the stability studies of these 
defects have shown that dynamical instabilities appear 
on the brane when there are more than one non-compact 
extra-dimensions [H^ . l53l l54| , whereas this is not the case 
for a five-dimensional bulk |55| provided that the usual 
fine-tuning between model parameters is fixed |5t| . 

The next step is now to derive observational con- 
sequences of braneworld cosmological models, e.g. the 
anisotropies of the cosmic microwave background (CMB). 
To that end, a lot of work has recently been in- 
vested to derive gauge invariant perturbation theory in 
braneworlds with one co-dimension [5?l l58l I59I l6fl | . Ac- 
cording to the background cosm olog y, the perturbations 
can be derived for a brane at rest ]6dJ , or moving in a per- 
turbed anti-de Sitter spacetime EE EE 0, EE EH • 
Whatever the approach chosen, the perturbation equa- 
tions are quite cumbersome and it is difficult to extract 
interesting physical consequences analytically. Also the 
numerical treatment is much harder than in usual four 
dimensional perturbation theory, since it involves partial 
differential equations. 

Nevertheless, it is useful to derive some simple physi- 
cal consequences of perturbation theory for braneworlds 
before performing intensive numerical studies. This has 
been done for tensor perturbations on the brane in a 
very phenomenological way in Ref. [(57j or on a more fun- 
damental level in Ref. 681. Tensor modes in the bulk 
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which induce scalar perturbations on the brane have been 
studied in Ref. |69l | and let to important constraints for 
braneworlds. 

In this article we consider a braneworld in a five dimen- 
sional bulk where cosmology is induced by the motion of 
a "3-brane" in AdS 5 . The bulk perturbation equations 
are considered without bulk sources and describe grav- 
ity waves in the bulk. The present work concentrates on 
the part of these gravity waves which results in vector 
perturbations on the brane. 

For the sake of clarity, we first recall how cosmology 
on the brane can be obtained via the junction conditions, 
particul arly emphasizing how Zi symmetry is imple- 
mented jlalllllillillillllllil- After re-deriving the 
bulk perturbation equations for the vector components in 
terms of gauge invariant variables [34 . l62l l63l l64l l65| , we 
analytically find the most general solutions for arbitrary 
initial conditions. The time evolution of the induced vec- 
tor perturbations on the brane is then derived by means 
of the perturbed junction conditions. The main result of 
the paper is that vector perturbations in the bulk generi- 
cally give rise to vector perturbations on the brane which 
grow either as a power law or even exponentially with re- 
spect to conformal time. This behavior differs essentially 
from the usual decay of vector modes obtained in stan- 
dard four-dimensional cosmology, and may lead, at least, 
to observable effects of extra-dimensions in the CMB. 

In the next section, the cosmological braneworld model 
obtained by the moving brane in an anti-de Sitter bulk 
is briefly recalled. In Sect. IIIII we set up the vector 
perturbation equations and solve them in the bulk. In 
Sect. HVl the induced perturbations on the brane are de- 
rived and compared to the usual four-dimensional ones, 
while Sect. E] deals with the consequences of these new 
results on CMB anisotropies. The resulting new con- 
straints for viable braneworlds are discussed in the con- 
clusion. 



II. BACKGROUND 



Another coordinate system for anti-de Sitter space 
which is often used in braneworld models is defined by 
r 2 /L 2 = e~ 2e / L so that 



ds 2 = g AB dx A dx B = e- 2e ' L (-dt 2 + <), ,d.r' d.r' ) + dg 2 



(4) 



A. Embedding and motion of the brane 



The position of the brane in the AdS 5 bulk is given by 



(5) 



where X M are embedding functions depending on the in- 
ternal brane coordinates (/i = 0, • • • , 3). Making use 
of reparametrization invariance on the brane, we choose 
x % = X % = y l . The other embedding functions are writ- 
ten 



X°=t h (r), X 4 = r h (r), 



(6) 



where r = y° denotes cosmic time on the brane. Since 
we want to describe a homogeneous and isotropic brane, 
X° as well as X 4 are required to be independent of the 
spatial coordinates y % . The four tangent vectors to the 
brane are given by 



Ma _ 9X M 

^ ~ dy» 



(7) 



and the unit space-like normal 1-form n M is defined (up 
to a sign) by the orthogonality and normalization condi- 
tions 



n M e" = 0, g AB n A n B = 1. 



(8) 



Adopting the sign convention so that n points in the di- 
rection in which the brane is moving (growing for an 
expanding universe), one finds using 



As mentioned in the introduction, our universe is con- 
sidered to be a 3-brane embedded in five dimensional 
anti-de Sitter spacetime with flat three dimensional spa- 
tial sections 



ds 2 = g AB dx A dx B = 



L 2 



-dt 2 



dij dx l dx°) 



L 5 



dr 2 . 



(1) 

The capital Latin indices A, B run from to 4 and i,j 
from 1 to 3. Anti-de Sitter spacetime is a solution of Ein- 
stein's equations with a negative cosmological constant A 



e° T =t b , e 4 T =r h , e\ (9) 

the components of the normal 

n = -r' b , 7i 4 = *b) n i = 0- (10) 

The other components are vanishing, and the dot denotes 
differentiation with respect to the brane time r. 

This embedding ensures that the induced metric on 
the brane describes a spatially flat homogeneous and 
isotropic universe, 



G AB + Ag AB = 0, 
provided that the curvature radius L satisfies 

*-4 



(2) 



(3) 



ds 2 = q^dy^dy" = -dr 2 + a 2 (r)%dyMz/, (11) 

where a(r) is the usual scale factor, and q^ v is the pull- 
back of the bulk metric onto the brane 



(12) 
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(see e.g. |jj 
defined by 



71|). The first fundamental form q AB is now 



(13) 



i.e. the push-forward of the inverse of the induced metric 
tensor |7(J, |72( • One can also define an orthogonal pro- 
jector onto the brane which can be expressed in terms of 
the normal 1-form 



n A n B — 9A1 



(14) 



in the case of only one codimension. 

Upon inserting the equations (JJJ, 1)1 0|l and (|13J1 into 
the above equation, one finds a parametric form for the 
brane trajectory [32. IM l6l l65| 



Th{r) = a(r)L, 

t h (r) = -y/T+ZMP, 



(15) 



where H = a/a denotes the Hubble parameter on the 
brane. Alternatively, comparing expression O with the 
Fricdmann metric Hll() yields the same result. 

Therefore, the unperturbed motion induces a cosmo- 
logical expansion on the 3-brane if is growing with 



with M 5 and G 5 , the five dimensional (fundamental) 
Planck mass and Newton constant. The superscripts ">" 
and "<" stand for the bulk sides with r > rb and r < rv 
As already noticed, the brane normal vector n M points 
into the direction of increasing r [see Eq. I|10fl ]. Eq. I|17|) 
is usually referred to as second junction condition. The 
first junction condition simply states that the first fun- 
damental form (|13fl is continuous across the brane. 

In general, there is a force acting on the brane which is 
due to its curvature in the higher dimensional geometry. 
It is given by the contraction of the brane energy momen- 
tum tensor with the average of the extrinsic curvature on 
both sides of the brane 28] 



S» v (K>, + K<„) =2/ = 0. 



}1V 



(19) 



This force /, normal to the brane, is exerted b y th e asym- 
metry of the bulk with respect to the brane [23, |Z£J ■ In 
this paper, we consider only the case in which the bulk 
is Zi symmetric across the brane, so that / = 0. In 
this case the motion of the brane is caused by the stress 
energy tensor of the brane itself which is exactly the cos- 
mological situation we have in mind. 

From Eqs. {TOl, dJ, {EJ) an d 03, noting that the 
extrinsic curvature can be expressed purely in terms of 
the internal brane coordinates 64, 631, one has 



B. Extrinsic curvature and unperturbed junction 
conditions 

The Lanczos-Sen-Darmois-Israel junction conditions 1 
govern the evolution of the brane. They relate the jump 
of the extrinsic curvature across the brane to its surface 
energy-momentum content jlj [Z~3- EEl EH ■ The extrinsic 
curvature tensor projected on the brane can be expressed 
by 



(16) 



Here V denotes the covariant derivative with respect to 
the bulk metric, and C n is the five dimensional Lie- 
derivative in the direction of the unit normal on the 
brane. With the sign choice in Eq. (|lf>|> . the junction 
conditions read [77j 

K^ v — K^ v = k 5 ( S^v — —Sq^u j = k 5 5*/j^, (17) 



K ihV = -\- [g AB (e*d„n B + e^d^n 3 ) + e*e*n c g AByC ] . 

A short computation shows that the non- vanishing com- 
ponents of the extrinsic curvature are 



1 + L 2 H 2 +L 2 H 



LVl + L 2 H 2 
Ku = - — Vl + L 2 H 2 Si 



(21) 



It is clear, that the extrinsic curvature evaluated at some 
brane position does not jump if the presence of the 
brane does not modify anti-de Sitter space. Like in the 
Randall-Sundrum model [l6| . in order to accommodate 
cosmology, the bulk space-time structure is modified by 
gluing the mirror symmetric of anti-de Sitter space on 
one side of the brane onto the other |3J|. There are two 
possibilities: one can keep the "r > rb" side and replace 
the "r < rb" side to get 



where is the energy momentum tensor on the brane 
with trace S, and 



k 2 = 6tt 2 G k 



1 

M?' 



(18) 



K > — K K < — — K 



(22) 



where is given by Eq. (|2ip. Conversely, keeping the 
r < r b side leads to 



K > = —K 
- Ii /if - fl /«f, 



K K = K 
zl /if J1 /i^- 



(23) 



1 In the following, they will be simply referred to as 
conditions" . 



'junction 



Note that both cases verify the force equation Q19[l . From 
the time and space components of the junction conditions 
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(|TT|l one obtains, respectively 



± 



1 + L 2 H 



= -K 2 (P + 0) - — K (P+T), 

L 2 H 2 2 sV ' 6 5 



± 



v/i + l 2 # 2 
I 



G 



«2 G° + T )- 



(24) 
(25) 



Here the brane stress tensor is assumed to be that of a 
cosmological fluid plus a pure tension T, i.e. 



Spu = (P + p) u^u v + Pq^ - Tci 



(26) 



p and P being the usual energy density and pressure on 
the brane, and the comoving four-velocity. The "±" 
signs in Eqs. (|21Jl and (|2*S|l are obtained by keeping, re- 
spectively, the r > rb, or r < ?*b, side of the bulk. In 
order to allow for a positive total brane energy density, 
p+T, we have to keep the r < n> side and glue it symmet- 
rically on the r > rb one 2 . In the trivial non-expanding 
case this construction reproduces the Randall-Sundrum 
solution with warp factor exp(— \g\/L), for — oo < g < oo. 
In our coordinates, we just have < r < rb on either side 
of the brane, and the bulk is now described by two copies 
of the "bulk behind the brane" . Even if r only takes val- 
ues inside a finite interval, and even though the volume 
of the extra dimension, 



V 



is finite, the bulk is non-compact and its spectrum of per- 
turbation modes has no gap (like in the Randall-Sundrum 
model). 

From Eqs. (|24|l and l|25|) . one can check that energy 
conservation on the brane is verified 

p + 3H (P + p) = 0. (28) 

Solving Eq. (|25|1 for the Hubble parameter yields 



H 2 



-5— p 1 
18 H 1 



P) + ^ T 2 

2TJ 36 



1 

L 2 ' 



(29) 



At "low energies", \p/T\ <C 1, the usual Friedmann equa- 
tion is recovered provided the fine-tuning condition 



36 



1 

T 2 '' 



(30) 



is satisfied. The four dimensional Newton constant is 
then given by 



K T 



(31) 



requiring a positive tension in order to get a positive 
effective four-dimensional Newton constant. Note also 
that low energy simply means r 2 ~ H~ 2 ^> L 2 . In the 
Fricdmann-Lemaitre era, the solution of Eq. I|29|l reads 



h~hA — 



-3(l+™)/2 



3 , o / a 



-3(l+w) 



(32) 



for a cosmological equation of state P — wp with con- 
stant w. The parameters H and a refer, respectively, 
to the Hubble parameter and the scale factor today. For 
the matter era we have w — 0, and during the radiation 
era w = 1/3. 



III. GAUGE INVARIANT PERTURBATION 
EQUATIONS IN THE BULK 

A general perturbation in the bulk can be decomposed 
into "3-scalar" , "3- vector" and "3-tensor" parts which 
are irreducible components under the group of isometries 
(of the unperturbed space time) 5*0(3) x E 3 , the group 
of three dimensional rotations and translations. In this 
paper we restrict ourselves to 3-vector perturbations 3 and 
consider an "empty bulk" , i.e. the case where there are 
no sources in the bulk except for a negative cosmological 
constant. Therefore we consider only bulk gravity waves 
since they are the only modes present when the energy 
momentum tensor is not perturbed. It is well known 
(see e.g. Ref. 0) 

that gravity waves in 4 + 1 dimensions 
have five degrees of freedom which can be decomposed 
with respect to their spin in 3 + 1 dimensions into a spin 
2 field, the ordinary graviton, a spin 1 field, often called 
the gravi-photon and into a spin field, the gravi-scalar. 
In this work we study the evolution of the gravi-photon 
in the background described in the previous section. 

After setting up our notations, we define gauge invari- 
ant vector perturbation variables in the bulk and write 
down the corresponding Einstein equations. Analytic so- 
lutions for all vector modes in the bulk are then derived. 



A. Bulk perturbation variables 

Considering only vector perturbations in the bulk, the 
five dimensional perturbed metric can be parameterized 



2 Note that we obtain the same result as in Ref. f>^, i.e. a positive 
brane tension for an expanding universe is obtained by keeping 
the anti-de Sitter side which is behind an expanding brane, ac- 
cording to its motion. 



3 For sake of clarity, the prefix "3-" will be dropped in what follows, 
and the term "vector" will be always applied here for spin 1 with 
respect to the surfaces of constant t and r. 
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dS 2 = -jz<it 2 + — (% + ViEj + VjSi) dx'dx^ 

T2 2 

+ — dr 2 + dMa; 4 + 2C i dx i dr, 

r 2 ^ 



(33) 



where Vj denotes the connection in the three dimensional 
subspace of constant t and constant r. Assuming this 
space to be flat we simply have Vj = di. The quantities 
E l . B l , and C % are divergenceless vectors i.e. c^-E^ = 
d.B 1 = d t C l = 0. 

As long as we want to solve for the vector perturbations 
in the bulk only, the presence of the brane is not yet 
relevant. Later it will appear as a boundary condition 
for the bulk perturbations via the junction conditions 
that will be discussed in Sect. II V Cl 

Under a linearized vector type coordinate transforma- 
tion in the bulk, x M — > x M + e M , with e M = (0, £», 0), the 
perturbation variables defined above transform as 



Ei 
Ci 



L 2 
L 2 



Ei + —Si, 



B t + —d t Si 
Ci + d r £i 



(34) 



As expected for three divergenceless vector variables 
and one divergenceless vector type gauge transformation, 
there remain four degrees of freedom which can be chosen 
as the two gauge invariant vectors 



Hi — Bi 
Sj = Ci 



d t E h 

,.2 



jd r Ei 



(35) 
(36) 



Note that in the gauge Ei = these gauge-invariant vari- 
ables simply become Bi and d respectively. 



B. Bulk perturbation equations and solutions 

A somewhat cumbersome derivation of the Einstein 
tensor from (|33|l to first order in the perturbations leads 
via Eqs. J3J) to the following vector perturbation equa- 
tions, 



0, (37) 



^d?Z + b^d r Z-L 2 d 2 Z + L 2 AZ = 0, (38) 



L 2 r \L 3 



L 2 \L 3 



+L 2 A[ ) = 0, (39) 



where A denotes the spatial Laplacian, i.e. 

A = piQidj, (40) 



and the spatial index on the perturbations has been omit- 
ted. One can check that these equations are consistent, 
e.g. with the master function approach of Ref. ■ 

A complete set of solutions for these equations can eas- 
ily be found by Fourier transforming with respect to x\ 
and making the separation ansatz: 

H(i,r,k) = 3 T (t,k)3 R (r,k), (41) 
E(t,r,k) = E T (t,k)E R (r,k). (42) 

The most general solution is simply a linear combination 
of such elementary modes. Eq. 1|38|) now leads to two 
ordinary differential equations for Z T and S R , 



I U r i — i R o U r ^ R 



^ + k 2 = ±n 2 , 



(43) 
(44) 



where k is the spatial wave number, and ±f2 2 the separa- 
tion constant having the dimension of an inverse length 
squared. We choose ft 2 > and indicate the two signs 
in Eqs. (|4*3*|) and i|4"4"fl by "±". The frequency O repre- 
sents the rate of c hange of 5 R at r ~ L, while the rate of 
change of S T is ^/|0 2 =p k 2 \. Eq. (|4*3*)) is a Bessel differen- 
tial equation of order two for the "— " sign and a modified 
Bessel equation of order two for the "+" sign [7SI |. while 
Eq. (|44|l exhibits, respectively, oscillatory or exponential 
evolution with respect to the bulk time. From Eq. QS) , 
similar equations are derived for E T (t, k) and E R (r, k). 
This time, the radial function is given by Bessel func- 
tions of order one. The constraint equation (|37(l ensures 
that the separation constant ±J7 2 is the same for both 
vectors and it also determines the proportionality factor 
of the amplitudes. The general solution of Eqs. I|37l) to 
ijjSH is a superposition of modes (fi, k) which are given 
by 



c oc < 



k 2 



^jj^n) e±^ 



(45) 



L -YS L -*\ e±^ 



Tv/T^fcW^ii — n e 



l 2 ^ ( £ 2 „\ ± t ^/n 2 - k 2 



E oc { 



±yji + k 2 /ti 2 —j 1 [ —n I /. 



±ity/n 2 + k 2 



(46) 
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Here K p and I p are the modified Bessel functions of 
order p while J p and Y p are the ordinary ones. The ± 
signs in Eqs. (|45(l and (|46(l come from the two indepen- 
dent solutions stemming from the second order differ- 
ential equation l|44|l . which is independent of the sign of 
the separation constant. In general, each of these modes 4 
can be multiplied by a proportionality coefficient which 
depends on the wave vector k and the separation con- 
stant ±f2 2 and which has to be the same for S and S. 
For f} 2 > k 2 the K- and I-modes can have an expo- 
nentially growing behaviour whereas for fl 2 < k 2 we use 
\/Q 2 — k 2 = i\/\n 2 — k 2 \ such that they become oscilla- 
tory. The J- and Y-modes are always oscillating. 

For a given perturbation mode to be physically accept- 
able one has to require that, at some initial time t\, the 
perturbations are small for all values < r < n>(fi), 
compared to the background. To check that, we use the 
limiting forms of the Bessel functions [j^. For large ar- 
gument, the ordinary Bessel functions behave as 



J p (x) 



Y p (x) 



I 2 ( 7T 7T\ 

\l^ COS { X -2 p -4)> 

I 2 / 7T 7T 

\ / — sin x — —p - 

V ttx V 2 l 4 



(47) 



while the modified Bessel functions grow or decrease ex- 
ponentially 



K p (x) 
l p (x) 



2x 
1 



(48) 



It is important to note at this point that, in this sense, 
also the J-modes and Y-modes are normalizable. One 
might view this integrability condition as a requirement 
to insure fmitcncss of the energy of these modes. This 
suggest that the J- and Y-modes could also be excited 
by some physical process. Indeed, from Eq. |0J, their 
divergence for r — > can be recast in terms of the g 
coordinate, with g — > oo. Expressed in terms of g, the 
integrability condition i|51[l ensures that the J- and Y- 
modes are well defined in the Dirac sense, and thus that 
only a superposition of them may represent physical per- 
turbations 5 . However, at small r, even if integrable, the 
perturbation amplitudes of the J- and Y-modes become 
large and it is not clear to us in which sense linear grav- 
itational perturbation theory can remain valid for these 
modes. 

During a bulk inflationary phase which leads to the 
production of 4 + 1 dimensional gravity waves, one may 
expect the K-modes to be generated, and for a given in- 
flationary model, it should be possible to calculate the 
spectrum of fluctuations, |^4(k, f2)| 2 . Clearly these fluc- 
tuations dominate the others and are extremely danger- 
ous since exponentially growing modes could rapidly lead 
to large fluctuations and destroy the homogeneity and 
isotropy of the universe. 

Nevertheless, in order to conclude about the viability 
of these modes, we first have to derive the evolution of 
their observational counterparts on the brane, especially 
in terms of cosmic time. This is done in the next section 
through the perturbed junction conditions. 



Therefore, in Eqs. I|45l) and (|46|l all modes, except for 
the K-mode, diverge as r — ► 0. Hence the only regular 
modes are 



S = A(k, 0)- 



£ = ±A(k,CL)\ 1 



k 2 i . (L\..\ ±fv /^r 



(49) 



-Ki 



(50) 

where the amplitude A(k, f2) is determined by the initial 
conditions and carries an implicit spatial index. For small 
wave numbers, k 2 < fl 2 , we are of course only interested 
in the growing mode but for large wave numbers both 
modes are comparable and oscillating in time. It is easy 
to see that the K-mode is also normalizable in the sense 
that 

2 



/ VWm 2 dr (x / - 
Jo Jo r 



1 1 

K,i — n 



ff ||E| 2 droc 



r 



7-2 
r 



dr < oo, 

(51) 

dr < oo. 



IV. GAUGE INVARIANT PERTURBATION 
EQUATIONS ON THE BRANE 

A. Brane perturbation variables 

Since we are interested in vector perturbations on the 
brane induced by those in the bulk, we parameterize the 
induced perturbed metric as 



dS 2 = q^dy^dy" 



-dr + 2abidrdy l 



(52) 



a 2 (5ij 



Vie, + Vj-eOdyM^, 



where e l and b % are divergenceless vectors. The junction 
conditions which relate the bulk perturbation variables to 
the perturbations of the brane can be written in terms of 
gauge invariant variables. Under an infinitesimal trans- 



formation yV 



where ^ = (0, a 2 ^), we have 



aii 



(53) 



4 In the following, they will be labelled by the kind of Bessel func- 
tion they involve, e.g. "K— mode" , "I— mode" . . . 



This is of course not the case for the I-modes. 
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Here the dot is the derivative with respect to the brane 
time r. Hence the gauge invariant vector perturbation 

is m 

(Ti = bi- adi. (54) 

It fully describes the vector metric perturbations on the 
brane. 

The brane energy momentum tensor S^ v given in Eq. 
IPS} has also to be perturbed. As we shall see, the junc- 
tion conditions (together with Z 2 -symmetry) do require a 
perturbed energy momentum tensor on the brane. Since 
we only consider vector perturbations 5p = SP = 0. 
However, the four-velocity of the perfect fluid does con- 
tain a vector perturbation = + du^ , with 




(55) 



and where v l is divergenceless. Under y M — > y^ + 

Vi^Vi-a^i, (56) 

where t>, = SijV 1 . A gauge invariant perturbed velocity 
can therefore be defined as 6 

& i =Vi + a&i. (57) 

In addition, the anisotropic stresses contain a vector com- 
ponent denoted m. Since the corresponding background 
quantity vanishes, this variable is g aug e invariant accord- 
ing to the Stewart-Walker lemma |8Cj . 

In summary, there are three gauge invariant brane per- 
turbation variables. We shall use the combinations 



Oi = h 

■di = Vi 

7T,;. 



(58) 



To apply the junction conditions we need to determine 
the perturbations of the reduced energy momentum ten- 
sor defined in Eq. I|17f) . In terms of our gauge invariant 
quantities they read 



SS TT = 0, 



(59) 



<VS' r , - -<• ( P+^p-^T j ,T, -U(P +,,)<>;. (<>()) 



5Sij = a 2 P {diir,j + djiTi) 



(61) 



B. Perturbed induced metric and extrinsic 
curvature 

We now express the perturbed induced metric, and the 
perturbed extrinsic curvature in terms of the bulk per- 
turbation variables |G5| . In principle there are two con- 
tribution to the brane perturbations: perturbations of 



Alternatively, one may use the fact that the covariant component 
Ui is already gauge invariant to find the gauge invariant quantity. 



the bulk geometry as well as perturbations of the brane 
position. A bulk perturbed quantity has then to be evalu- 
ated at the perturbed brane position [see Eq. jEH- Using 
reparametrization invariance on the brane [65j , the latter 
can be described by a single variable T, 

X M = X M + Tn M , (62) 

where all quantities are functions of the brane coordi- 
nates y^ . Since T is a scalar perturbation it does not 
play a role in our treatment, and we can consider only 
the perturbations Sg AB due to the perturbed bulk geom- 
etry evaluated at the unperturbed brane position. The 
induced metric perturbation is given by 

= <7m" ~ = e ti e ^9AB- (63) 
From Eqs. (|33|) . (|3l))l and (|p^|l one finds in the gauge 
Ei=0 

5q TT = 0, 

Sq Ti = oVl + L 2 H 2 Ei + aLHY^i, (64) 



= <*9i 



0. 



The time component vanishes as it is a pure scalar, and 
the purely spatial components can be set to zero without 
loss of generality by gauge fixing. 

In the same way, perturbing Eq. (|20|l . and making use 
of Eqs. Q, 10 in order to derive the perturbed normal 
vector, leads to (again we use the gauge Ei — 0) 

(65) 



8K T 







[l + L z H 2 )Ei 



5K lj = -aLH(diZj + 9,-Sj) 



(66) 



(67) 



where SK.. 



K„ 



i/iu — j-^fj.u - K^ V1 and all bulk quantities have 
to be evaluated at the brane position. In the derivation 
we have also used that on the brane <9 M = e^d A . 



C. Perturbed junction conditions and solutions 

The first junction condition requires the first funda- 
mental form q AB to be continuous across the brane. 
Therefore, the components of the induced metric ijol'l) 
are given by the explicit expressions (|64l) . This leads to 
the following relations 



&i Ei , 

hi = Vl + L 2 H 2 B. t + LHQ, 



(68) 



where the bulk quantities have to be taken at the brane 
position (ib, J"t>)- For Oi — hi — aet we use 



= a (t h d t Ei + r h d r Ei) 

= \J\ + L 2 H 2 d t Ei + ^LHdrE, 



(69) 
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Together with Eqs. I|35[l and (|36|l this gives 

o-j = Vl + L 2 H 2 E, + LHY, Z . (70) 

The equations corresponding to the second junction con- 
dition are obtained by perturbing Eq. l(T7|) (using K> = 
—K< y = —Kpv) and inserting the expressions (|66[) . 
()67|l for the perturbed extrinsic curvature tensor, with 
Eqs. I(60|) . (|61(l for the perturbed energy-momentum ten- 
sor on the brane. After some algebra one obtains, respec- 
tively, for the (Oi) components, and the (ij) components 
respectively 

7=== ? « (*< + #i) = « 2 9rS, - d&i, (71) 

K 2 5 aP7Ti = -LHEi - ^l + L 2 H 2 Y, u (72) 

where we have used the unperturbed junction conditions, 
Eqs. H24|) and i|25|) . and the fact that on the brane = 
e^d A . One can use Eq. (|70[) to eliminate the variable 
in Eq. 1)72[) in favor of Oi . This leads to 

K 2 s aLHPiT 2 + \fi + L 2 H 2 a z = 3*. (73) 

Hence, if by some mechanism, like e.g. bulk inflation, 
gravity waves are produced in the bulk, their vector parts 
Sj and Si will induce vector perturbations <Ji, vorticity 



Oi + $i , and anisotropic stresses ~ki on the brane accord- 
ing to Eqs. J23), jnj and These are necessary to 
satisfy the boundary conditions implied by Z2 symmetry. 
In general, Ui, Ui+di, and 7r^ come from both, bulk per- 
turbations and intrinsic perturbations on the brane, but 
here, we want to study only the perturbations induced 
from bulk gravity waves and do not consider additional 
sources on the brane. Then the brane variables are fully 
determined by Sj and S^. 

In the following we do not want to specify a particu- 
lar mechanism which generates Sj and Sj, and just as- 
sume they have been produced with some spectrum given 
by A(k, J7). In 4-dimensional standard cosmology it is 
well known that vector perturbations decay, and there- 
fore they are not considered. Note that in simple in- 
flationary models they are not even generated. Only if 
they are continuously re-generated like, e.g. in models 
with topological defects (see e.g. Ref. they affect 

CMB anisotropics. Here the situation is different since 
the modes considered are either exponentially growing 
or oscillating at least with respect to the bulk time. To 
calculate the CMB anisotropics from the vector pertur- 
bations induced by bulk gravity waves, the relevant quan- 
tities are a and 1? + <J given in terms of the bulk variables 
by Eqs. (|70|) and (|71|) . Inserting the solutions (|49|l and 
(j5U|l for the K-mode into J7UJ and (fTTf 1 yields 



a(t h ,k) = A(k,il) 

((7 + 1?) (*b,k) 



/ rrrl / LQ\ / k 2 1 / Lfl\ 



±t b yJW-k 2 



A(K^ VTTim "^K 1 (™) e±^ 
^ 2 2L 2 H a a 2 \ a J 



(74) 
(75) 



where again we have omitted the spatial index i on a, 
i? and A. Similar equations can be obtained for the J- 
and Y-modes by replacing, in Eqs. I)74|) and (|75|l . the 
modified Bessel function by the ordinary ones, plus the 
transformations: — k 2 — > k 2 , H — > —H and t — > it. 

These equations are still written in bulk time which 
is related to the conformal time r\ on the brane by 



dt h = Vl + L 2 H 2 dr}. 



(76) 



Therefore, at sufficiently late time L 2 H 2 -C 1 such that 
dtb — dr/. Note that L is the size of the extra-dimension 
which must be smaller than micrometers while H^ 1 is 
the Hubble scale which is larger than 10 5 light years at 
times later than recombination which are of interest for 
CMB anisotropics. 

As a result, the growing or oscillating behaviour in bulk 
time remains so in conformal time. Moreover there are 
additional time dependent terms in Eqs. (|74|l and l(75|l 
with respect to Eqs. 149|) and (|50|l due to the motion of 



r 



the brane. As can be seen from Eqs. 1(74(1 and 1)75(1 . the 

modes evolve quite differently for different values of their 
physical bulk wave number Q/a. In the limit f2/a <C 1/L 
and for £l 2 > k 2 , the growing K-modes behave like 



2A 

(nZ) 5 

A 



k 2 



(ML) 2 2a 2 H 



(77) 



where use has been made of L 2 H 2 -C 1, and of the lim- 
iting forms of Bessel function for small arguments [T^ 



1 



x^O z 



(78) 



In the same way, from Eq. I|48f) . the K-modes verifying 
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ft /a ^> 1/L reduce to 



CMB ANISOTROPIES 



A 



(my 

A 



k 2 



(ML) 2 2a 2 H 



,t i <>/. v ; ' 

2 V a J 



ClL/a 



1/2 



a 



1 _ 



They are exponentially damped compared to the former 
[see Eq. (|77|) ]. As a result, the main contribution of the 
K-mode vector perturbations comes from the modes with 
a physical wave number fi/a smaller than the energy 
scale 1/L associated with the extra-dimension. As the 
universe expands, a mode with fixed value remains rel- 
atively small as long as the exponents in Eq. I|79|l satisfy 

—L - rj^/n 2 -k 2 ~-(riL- r^/n 2 - k 2 ) > o. (so) 

a a \ I 

When this inequality is violated, for k <C this is soon 
after r ~ L, the mode starts growing exponentially. The 
time t ~ L also corresponds to the initial time at which 
the evolution of the universe starts to become Friedman- 
nian. 

In the same way, one can derive the behaviors of the 
J- and Y-modes on the brane for physical bulk wave 
numbers greater or smaller than the size of the extra- 
dimension. This time, the exponentially growing terms 
are replaced by oscillatory ones, and the ordinary Bessel 
functions are approximated by {see Eq. l(T7|) and |7q} 



3 p (x) 



z-»o r(p + 1) 



1 /xy 



1 /2 XP 

y p (x) ~ — r(p) 



(81) 



From Eqs. JSJ and J^U, the equivalents of Eqs. J73J 
and 1)75(1 for J- and Y-modes can be shown to oscillate 
always. From Eq. their amplitude is found to decay 
like a~ 3 ' 2 in the short wavelength limit Q/a^§> 1/L. In 
the long wavelength limit VI /a <C 1/L, the amplitude of 
the Y-mode stays constant whereas the J-mode decreases 



as a 



The vorticity is also found to oscillate in conformal 
time. This time, the amplitude of the long wavelength 
Y-modes always grows as a 3w+1 while for the J-modes 



it behaves like 



,3u;-l 



Finally, in the short wavelength 



limit, both Y and J vorticity modes grow like a 3w+1 ^ 2 . 

Whatever the kind of physical vector perturbation 
modes excited in the bulk, we have shown that there 
always exist bulk wave numbers f2 that give rise to grow- 
ing vector perturbations on the brane. Although the J- 
and Y-modes generate vorticity growing like a power law 
of the scale factor, they can be, in a first approximation, 
neglected compared to the K-modes which grows like an 
exponential of the conformal time. In the next section, 
we will derive constraints on the amplitude A(k, fl) of 
the K-modes by estimating the CMB anisotropics they 
induce. 



In order to determine the temperature fluctuations 
in the CMB due to vector perturbations on the brane, 
we have to calculate how a photon emitted on the last 
scattering surface travels through the perturbed geome- 
try H52[> . A receiver today therefore measures different 
microwave background temperatures T R (n l ) for incident 
photons coming from different directions n' '. In terms of 
conformal time the vector-type temperature fluctuations 
are given by 



= n l (a t +dA\l + 



dxi' 



? dA, 



(82) 



where A denotes the affine parameter along the photon 
trajectory and the prime is a derivative with respect to 
conformal time 77. The "R" and "E" index refer to the 
time of photon reception and emission, which in the fol- 
lowing correspond to the conformal times today rj and 
at recombination r) E . For the second equality we have 
used 



d(7j 

dA 



= a, — n J 



dxi' 



(83) 



where — n % is the direction of the photon momentum. 
We have also neglected the contribution from the upper 
boundary, "R" , in the first term since it simply gives rise 
to a dipole term. The first term in Eq. (|82|l is a Dopplcr 
shift, and the second is known as integrated Sachs- Wolfe 
effect. To determine the angular CMB perturbation spec- 
trum Ci, we apply the total angular momentum formal- 
ism developed by Hu and White |82|. According to this, 
a vector perturbation v is decomposed as 



where 



e- = ^(e( 1 )±ie( 2 ) 
y/2 



(84) 



(85) 



and e^ 1,2 ) are defined so that (e^, e^ 2 \ k = k/fc) form 
a righthanded orthonormal system. Using this decom- 
position for $i and cr^, one obtains the angular CMB 
perturbation spectrum Ce via 



Ct = -M + l) k 2 (\A e (k)\ 2 )dk (86) 
I" Jo 



where 



A,(fc) = - l ?+(r 7E ,fc) 



j e (kr) n - kr) E ) 
kr] — krj 



(87) 



dry. 
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In Eq. (|87|l we have assumed that the process which gen- 
erates the fluctuations has no preferred handedness so 
that (|cr+| 2 ) = (|cr-| 2 ) as well as (|i9+| 2 ) = (\tf~\ 2 )- Omit- 
ting the "±" superscripts, we can take into account the 
negative helicity mode simply by a factor 2. 

As shown in the previous section, the main contribu- 
tion of the K-modes comes from those having long wave- 
lengths a /CI S> L, and k < CI. In the following, only these 
modes will be considered. Since they are growing expo- 
nentially in r), the integrated Sachs- Wolfe contribution 
will dominate and we concentrate on it in what follows. 
A more rigorous justification is given in the appendix. 
Inserting the limiting form l|77|) for a in Eq. (|87|l gives 



A e (k) ~ 2A H 



k' 2 



x / ■ — — e 

/o x 



nix) p^zj 



dx, 



where a simple power law ansatz has been chosen for the 
primordial amplitude 



^(^(k,^)! 2 ) = A (ci)ci 2 L 2 k n . 



(89) 



A dimensionless wave number k, and conformal time rj, 
have also been introduced as 



(90) 



in order to measure their physical counterparts in unit of 
the bulk wavelengths. The condition k < CI becomes now 
k < 1. The integration over rj in the integrated Sachs- 
Wolfe term is now transformed into an integration over 
the dimensionless variable x defined by 



x = k (r) - T}) = k (r? - rj) 



(91) 



In the following, we will use the approximation x E = 

k iv ~Ve) - kr) - 

By observing the CMB, one may naturally expect that 
the perturbations with physical wavelength greater than 
the horizon size today have almost no effect. In terms of 
our parameters, this means that the main contribution 
in the Cg comes from the modes verifying Cl/a > H , 
hence fj > 1. 

In the appendix, we derive a crude approximation for 
the angular power spectrum induced by the exponentially 
growing K-modes, in a range a little more constrained 
than the one previously motivated, namely 



ci 

<H n < — < 



L- 1 
1 + z ' 



(92) 



where z E is the redshift at photon emission which is taken 
to coincide with recombination, z E ~ 10 3 . In order to 
simplify the calculation, we do not want the transition 
between the damped K-modes (CI/ a > and the ex- 
ponentially growing ones (CI /a < L~ r ) to occur between 



the last scattering surface and today. This requirement 
leads to the upper limit of Eq. I|92l) . Moreover, in order 
to derive the Ci, we perform an expansion with respect 
to a parameter £/r} assumed small, and £ max refers to 
the multipole at which this approximation breaks down. 
The lower limit in Eq. (|92|) comes from this approxima- 



tion. Using the values L ~ 10 ''mm, H 
g max ~ 10 3 , and z E ~ 10 3 , one finds 
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29 , 



10 26 mm 1 < — < 1mm 1 . 



(93) 



The corresponding allowed range for the parameter rj 
becomes [see Eq. (IA.11I) ] 



10 3 < rJ Q < 10 29 . 



(94) 



Clearly the detailed peak structure on the CMB 
anisotropy spectrum would have been different if we had 
taken into account the oscillatory parts (k > ft) of the 
K-modes, as well as the Y- and J-modes, but here we 
are only interested in estimating an order of magnitude 
bound. As detailed in the appendix, for a scale invariant 
initial spectrum, i.e. n — —3/2, we obtain 



(95) 



From current observations of the CMB anisotropies, the 
left hand side of this expression is about 10~ 10 , and for 
£ ~ 10, one gets 



AM 



< 



[«7 D -51n(ij Q )] 



(96) 



From Eq. (|A.11|I and l|94|l . one find that the primordial 
amplitude of these modes must satisfy 

A (Cl) < e~ w \ for VL/a ~ lO^mm" 1 (97) 

and, more dramatically, 

A {fl) < e" 10 ' 9 , for n/a a ~ lmm" 1 . (98) 

As expected, the perturbations with wavelength closer 
to the horizon today (smaller values of fl) are less con- 
strained than smaller wavelengths [see Eq. (|97|l ]. More- 
over, one may expect that the bound l|98|l is no longer 
valid for Q/a > L ~ 1 / ( 1 + z E ) since the modes in Eq. ifFi^) 
start to contribute. However, the present results concern 
more than 20 orders of magnitude for the physical bulk 
wave numbers CI /a , and show that the exhibited modes 
are very dangerous for braneworlds. 

It seems that the only way to avoid these constraints 
is to find a physical mechanism forbidding any excitation 
of these modes. 
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VI. CONCLUSION 

In this paper we have shown that vector perturbations 
in the bulk generically lead to growing vector perturba- 
tions on the brane in the Friedmann-Lemaitre era. This 
behaviour radically differs from the usual one in four- 
dimensional cosmology, where vector modes decay like 
a~ 2 whatever the initial conditions. 

Among the growing modes, we have identified modes 
which are perfectly normalizable and lead to exponen- 
tially growing vector perturbations on the brane with re- 
spect to conformal time. By means of a rough estimate 
of the CMB anisotropies induced by these perturbations, 
we have found that they are severely incompatible with a 
homogeneous and isotropic universe by lighting up a fire 
in the microwave sky, unless their primordial amplitude 
is extremely small. 

No particular mechanism for the generation of these 
modes has been specified. However, we may expect that 
bulk inflation leads to gravitational waves in the bulk 
which do generically contain them. Even if they are not 
generated directly, they should be induced in the bulk 
by second order effects. Usually, these effects are too 
small to have any physical consequences, but here they 
would largely suffice due to their exponential growth [see 
Eqs. (|97|l and (J5SJ]. This second order induction seems 
very difficult to prevent in the models discussed here. 

It is interesting to note that this result is also linked to 
the presence of a non-compact extra-dimension which al- 
lows a continuum in the bulk frequency modes f2. A 
closer examination of Eq. (|44|l shows that the mode 
fl = 0, admits only J- and Y-mode behaviours. In a 
compact space, provided the first quantized value of f2 
is sufficiently large, one could expect the exponentially 
growing K-modes to be never excited by physical pro- 



cesses. Another more speculative way to dispose of them 
could be to study their causal structure: By setting k = 
in Eq. I|44|) . one might be tempted to view — f2 2 as a neg- 
ative mass squared, suggesting that these K-modes could 
physically represent bulk tachyons. This would motivate 
their causal decoupling from the subsequent evolution of 
the brane. In a more basic theory, which goes beyond 
our classical relativistic approach, these modes may thus 
not be present at all. However, this interpretation breaks 
down for modes with k > f2 which can again influence 
the brane light cone. 

It is important to retain that even if the K-modes can 
be eliminated in some way, the growing behavior of the 
Y- and J-modes remains. Although their power law 
growth is not as critical as the exponential growth of 
the K-modes, they should have significant effects on the 
CMB anisotropies. Indeed, they lead to amplified oscil- 
lating vector perturbations which are entirely absent in 
four-dimensional cosmology. This will be the object of a 
future study [83|. 

We therefore conclude that, if no physical mechanism 
forbids the generation of the discussed vector modes 
with time dependence cx exp(?7vTJ 2 — k 2 ), anti-de Sitter 
braneworlds, with non-compact extra-dimension cannot 
reasonably lead to a homogeneous and isotropic expand- 
ing universe. 
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APPENDIX: CMB ANGULAR POWER SPECTRUM 



In this appendix we first present a crude and then a more sophisticated approximation for the C^-power spectrum 
from the exponentially growing K-modes. As we shall see, at moderate values of i ~ 10-50, both lead to roughly the 
same bounds for the amplitudes which are also presented in the text. 



1. Crude approximation 

Here we start from Eq. (JSSJ. In the integral 

: MEle-^VV^ds, (A.l) 
x 

we replace jt, by its asymptotic expansion for small £, 

^ x) ^&m0m- (A - 2) 

This is a good approximation if either x E ~ krjo < 1/2 or x E (l/k 2 — l) 1 / 2 > 1. The integral of x then gives 

(|A^)| 2 )-^^(i-l) (A.3) 
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Integrating over k, we must take into account that our approximation is only valid for k < fc max = (£! 2 — r] 2 ) 1 / 2 . 
Since we integrate a positive quantity we certainly obtain a lower bound by integrating it only until k max . To simplify 

the integral we also make the substitution y = yl — k 2 . With this, and inserting our result (|A.3|I in Eq. J5BJ), we 
obtain 

^ £ Sl<^ n+3 f (! " y 2 ) n+l - 1,2 y"- 2l e 2 ^ y Ay. (A.4) 

For £ > 2, y 3 ~ 2f > 1 on the entire range of integration. Hence we have 

PCt Z fUo fi2 " +3 J\l - v 2 )"^" 1 / Wdy. (A.5) 
This integral can be expressed in terms of modified Struve functions |78| . In the interesting range, rj g 3> 1 we have 

\l - y r +i - 1/2 e^y d y ~ e 2 % r( ^+J 2 /2) . (A.6) 



Inserting this result in Eq 1|A.5|) we then finally obtain 



£ 2 C e > V2^A 2 n 2n+3 e 2 % 1±£— [±] = V2^A 2 e 2 ^o 111— [±\ , (A.7) 



2 2i In i 2 21 \fj 

where we have used Stirling's formula for T(£ + n + 1/2) and set n = —3/2 after the equality sign. 

In the next section we use a somewhat more sophisticated method which allows us to calculate also the Doppler 
contribution to the Ct's. For the ISW effect this method gives 



for n = —3/2. Until £ ~ 15 the two approximations are in reasonable agreement and lead to the same prohibitive 
bounds for A (fi). For £ > 15, Eq. (|A.8(1 becomes more stringent. 



2. Sophisticated approximation 

In Eq. (|88|l we have only considered the dominant contribution coming from the integrated Sachs- Wolfe effect. The 
general expression is obtained by inserting the solutions (|77|) for a and d in Eq. I|87|> . 



(A.9) 

To derive the first term we have used Eq. (|32|l in the matter era. The parameter 

P=6(l + z E )(^) 2 , (A.10) 

reflects the change in behavior of the modes, redshifted by z E to the emission time, which are either outside or inside 
the horizon today. It is important to note that the parameter H a /Cl completely determines the effect of the bulk 
vector perturbations on the CMB, together with the primordial amplitude A . Indeed, solving Eq. (|29|l in terms of 
conformal time, and using Eqs. (|28|l and 1)32(1 . yields i] Q ~ 2/(a Q H n ) in the Friedmann-Lemaitre era. Thus 



77 ~ 2-^, 77 ~ £-a. (A.ll) 

H /e 1 + z e H K ' 



, fi/Q„ 1 n/a 

' H ' % ~ l + z E H 

We now replace the spherical Bessel functions ji in the integrated Sachs- Wolfe term (ISW) using the relation |78j 

3i{x) = t &J W /2(4 (A.12) 
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In the ISW term the upper integration limit can be taken to be infinity as the contribution from x E to infinity can 
be neglected provided X\J\jk 2 — 1 > 1. This restriction is equivalent to k 2 < 1 — l/rj o which is verified for almost all 
values of k up to one, given that Tj Q varies in the assumed range (|94|l . We remind that for the exponentially growing 
K-mode k < ft, and hence < k < 1. This allows for the exact solution [84| 

[ .-'VM*-^ = ^j^Lf (|. f + f !*■) , (A.i3) 

where F is the Gauss hypergeometric function. In regard to the subsequent integration over k we approximate F as 
follows. For small values of k, F is nearly constant with value 1, at k = 0. As k — > 1 the slope of F diverges and it 
cannot be Taylor expanded anymore. However, by means of the linear transformation formulas |78j |. F can be written 
as a combination of hypergeometric functions depending on 1 — k 2 

f{ 1 - e - + i-t+ 3 -p)= r (' + CK£> f (- £ + 

U'2 + ' + 2' „(l 3\„/* 1\ U'2 + '2' 



2 2 / \2 2 



; r 



(A. 14) 



+ r ^HH " + ~' 

These in turn can be expanded around 1 — k 2 = 0, and gives 

F ~ 2 £+1 / 2 (l-^Vl-fc 2 ) • (A.15) 



These two approximations intersect at k = yl — l/^ 2 . In this way, we can evaluate the mean value of F by integrating 
the two parts over the interval [0, 1]. Thus, the hypergeometric function is replaced by 

Furthermore, the Gamma functions in (|A. 13|> can be approximated using Stirling's formula |78j 

r(£) _ i 



r(£ + 3/2) £3/ 2 ' 
Putting everything together and squaring Eq. (|A.9|) we obtain 



(A.17) 



f / p\ 2 e 2 T E Vi-fe 2 

iA,(fc)i 2 = 2 ^w (i-y J ( _ __ )3 (v[%-y) 



( If) _- )3/2 W [*fo„ - *,)] ^ 

The C/s are then found by integrating over all k- modes 

C e = -til + l)fl 3 [ k 2 \A e (k)\ 2 dk 

n Jo (A.19) 

= AA 2 l{l + l)Q 2n + 3 (C { p + Cf + Cf ') , 

where the C^ 1 correspond to the three terms in Eq. (|A.18|I . In the following we keep only the zero order terms in 
V a /v E - From Eqs. (|A.18jl . (|A.19|) one finds 



{ P = ^ (* " i) 2 ^^ [ J ^V 2 (H)] 2 M (A.20) 



c, 

V Jo 
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First, notice that if the argument is larger or smaller than the index, the Bessel functions are well approximated by 
their asymptotic expansions (|47|l and l)8ip. respectively. Therefore, we split the k- integral into two integrals over the 
intervals [0, kg] and [kg, 1], in each of which the Bessel function is replaced by its limiting forms. The transition value 
kg is given by kg ~ ^/v - I n the integral from kg to 1, the sin 2 (A;7j o ) is then replaced by its mean value 1/2 which is 
justified if the multiplying function varies much slower in k than the sine. To carry out the integration we make the 
substitution y 2 = 1 — fc 2 , and in order to simplify the notation we define the integral 



Z(a,M)= / y(l - y 2 ye 2 Wdy 



(A.21) 



In this way we can write Eq. (|A.20fl in the form 



C, 



(i) 



' —4 



1(0, yg, n - 3/2) - =-T(0, y t , n - 1/2) + ^1(0, y t , n + 1/2) 



r 2 (£ + 3/2) 



_ x 2£+i r 

'0 



l(y e , 1,1 + n- 1/2) - iL J( w , 1, 1 + n + 1/2) + Jjl( w , l,£ + n + 3/2) 



(A.22) 



Since yg = y 1 — fc 2 is very close to one, and the integrand is continuous in the interval [0, 1], integrals of the form 
X(yt, 1, v) can be well approximated by the mean formula 



l{y h l,v)~y{l- y ye 



(i - Vi) 



2n 

e e 



2(i/+l) 



(A.23) 



y=w 2 

For the integrals of the type X(0, yg, v) we distinguish between three cases: 

Case a: v > — 1. This case corresponds to a spectral index n > 1/2 in the first integral in Eq. (|A.20J| . We write 
X(0, y^, v) — T(0, 1, i/) — I(yg, 1, f). The solution of the latter is given by Eq. (|A.23|I . whereas the former can be 
solved in terms of modified Bessel and Struve functions 841 



Since our derivation assumes 77 > I, the large argument limit applies and we have 



(A.24) 



W3/2(2r? E )+L !/+ 3 /2 (277 E ) 



277 

e e 



/7T77 F 



(A.25) 



independently of the index v. 
Case b: v = — 1. Since the above expressions, Eq. 1A.24II . diverge for v = —1, we approximate the integral by 



l{0,yt,v) ~ e 



2 % 



yil-y 2 )- 1 ^ 



e 2?7 E In 



(A.26) 



We have checked that the numerical solution of X(0,?/£,^) agrees well with the approximation, provided yg is 
close to 1. 



Case c: v < — 1. We use the same simplification as in Eq. i|A.26(l . and now the integral yields 



l(0,y e ,v) ~ e 



2 ^ 



in 



y(i-y 2 Ydy 



27) 

e e 



2(1/ + 1) V>7 



2(i/+l) 



(A.27) 



For the particular value 77 = —3/2, Eq. l|A.22Jl contains terms X(0, y i, —3) and T(0, yg,— 2) which can be evaluated 
according to (|A.27J| . as well as a term X(0, yg, —1) for which we use (|A.26J| . The remaining three integrals over the 
interval [yg, 1] are evaluated by (|A.23(I . The result is 
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The parameter (3 is a constant of order unity, within in our approximation it is [3 — 1 — In 2 ~ 0.3. 
The second term C^p in Eq. (|A.18J| reads 



C, 



(2) 



/2 {kr] o )dk, 



(A.29) 



where only the zero order terms in rj E /r] o has been kept. Using the limiting forms for the Bessel function for arguments 
smaller and larger than the transition value ke, yields 



C, 



(2) 



1 



3£ 7/2^3/2 r(£ + 3/2) J Q 



dfe 



(A.30) 



dfc. 



For consistency with the derivation of Cj , we have assumed that the main contribution comes from the first integral, 
while the second one is small due to the oscillating integrand. Since fc^ <C 1, we can use again the mean formula to 
evaluate the first integral, and by the Stirling formula for T(£ + 3/2), Eq. (|A.30|) becomes 



a 



(2) 



127r l/2£ll/2 



£ 



2n+£+2 



J3t 



(A.31) 



Since I < T) Q , the spectrum is damped at large i, while the other terms can lead to the appearance of a bump, 
depending on the value of rj and n. 

The last terms C) reads 



C, 



(3) 



36£ 7 



(A.32) 



Splitting this expression in two terms over 1 — fc 2 , and using the substitution y 2 = 1 — fc 2 yields 



C. 



(3) 



1 



m 7 



T 0,l,n- 



J 0,l,n- 



(A.33) 



where T is given by Eq. (|A.21I) with rj E — > r? n . As before, these two integrals can be expressed in terms of modified 
Bessel and Struve functions |84j . From Eq. I|A.24|) . taking their limiting forms at large argument, and expanding the 
r function by means of the Stirling formula gives 



a 



(3) 



72£ 7 (n + £ + 3/2) (n + 1 + 1/2) V 2 £ 15 ? 2 36 V V 



7r e 2?? o e 



n+i+1/2 



(A.34) 



Clearly, cj 3 "* dominates over the others since it involves exp(2ry ) while Cji and appear only with fractional 
power of this factor, namely exp(^ ) and exp(?y /z E ). This is due to the fact that we are concerned with incessantly 
growing perturbations leading to the predominance of the integrated Sachs- Wolfe effect. 

Inserting Eqs. (|A.28|) . (|A.31I) and l|A.34|) for the particular value n = -3/2 into Eq. (|A.19fl gives the final CMB 
angular power spectrum 



2tt tt °] 4tt 



12 7 r 1 /2^3/2 



6z E 
-l 



12z B 



In 



V 



2@l 



t 2 



1 - 



24z„ 



i 



tt e 2 "o e~ e 
2 £V2 "36~ \r] 



24z E 
e-i 



(A.35) 
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